1 



A Model for the 3 He(d, p) 4 He Reaction at Intermediate Energies*) 

H. Kamada**) >***)>****), W. G16ckle***\ H. Witala*****) , S. Gojuki****) 

***) Institut fur Theoretische Physik II, Ruhr- Universitdt Bochum, D-44^80 

Bochum, Germany 

****) D er p ar i men i j Physics, Faculty of Science and Technology, Science University 
"O ! of Tokyo, Noda, 278-8510 Chiba, Japan 

*****) l ns tit u t e f Physics, Jagellonian University, PL-30059 Cracow, Poland 

(Received February 8, 2008) 

D Polarization correlation coefficients have been measured at RIKEN for the 3 He(d,p) He 

reaction at intermediate energies. We propose a model for the (d, p) reaction mechanism 
' using the pd elastic scattering amplitude which is rigorously determined by a Faddeev cal- 

culation and using modern NN forces. Our theoretical predictions for deuteron polarization 
, observables A y , A yy , A xx and A xz at £^=140, 200 and 270 MeV agree qualitatively in shape 

^ ' with the experimental data for the reaction 3 He(d,p) 4 He. 

m : 
o 

O ' §!• Introduction 

<n : 

\ The measurement ^ of the 3 He (d, p) 4 -ffe reaction at RIKEN aimed at an in- 

vestigation of the high momentum components of the deuteron wave function and 
the d-state admixture linked to them. High precision data resulted for the polariza- 
^T* ! tion observables A y , A yy , A xx ,C y:V and C XjX . Out of them the linear combination 

q ■ Cm = 1 + \{A yy + A xx ) + \{C y ^y + C X)X ) has been formed x ). The Dubna and Sat- 

urne groups also obtained the polarization correlation coefficient Cm built in this 
case from the measurements of T20 and kq in d + p backward scattering > and from 
.^J I the inclusive deuteron breakup process 3 ). The polarization correlation coefficient C11 

^ ■ at forward angles of the outgoing proton is directly related to the ratio of deuteron 

wavefunction components if one uses the plane wave impulse approximation (PWIA) 

C|| PWIA) 4^2^^) + w 2 (k pn ) 1 ^ 

Here u and w are the S-, D-wave components of the deuteron wavefunction, and k pn 
a kinematically fixed relative momentum of the pn pair. These PWIA calculations 
are very poor in relation to the data x ) . This is shown in table I for Ci . There we also 
exhibit the different D-state probabilities for the modern realistic NN potentials, CD- 
Bonn 4 ), AV18 5) and Nijmegen 1,11 and 93 6) . Clearly one needs a better calculation 
for the analysis of the 3 .£/e((f,p) 4 He reaction. 
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Letters 



A theoretical analysis has been reported by the SUT group 7 ) based on a 3 He-n- 
p and d-d-p three-cluster model. However, the evaluations performed up to now in 
this model lead only to a tiny deviation from the PWIA calculations just mentioned. 
Recently, the Hosei group 8 ) analyzed T20 and kq with the 3 He-n-p cluster model 
by an analogy between 3 He and the proton (T=l/2, S=l/2). They conclude that 
PWIA describes the global features of the experimental data. 

In this letter we would like to introduce again a 3N model, which when evaluated 
correctly leads to a great similarity of various polarization observables to the ones 
found in the reaction 3 .£/e(<i,p) 4 He. 



Table I. The polarization correlation coefficient Cii according to Eq. (IT) in a simple PWIA model 
for different NN potentials and from an experimental value for the reaction 3 f?e(d,p) 4 He. We 
also show the corresponding deuteron D-state probabilities. 



Potential 


C\\{PWIA) 


D-state Probability (%) 


CD-Bonn 4 ' 


0.645 


4.86 


AV18 5) 


0.722 


5.78 


Nijmegen 93 6) 


0.710 


5.76 


Nijmegen I 6 ' 


0.712 


5.68 


Nijmegen II 6 ' 


0.726 


5.65 


exp. 4 ' 


0.223T 0.044(statistical)± 0.037(systematic) 





§2. Model 

For the 3 He(<i,p) 4 He reaction we assume a model which is based on a three-body 
reaction process. This is shown in Fig. 1. The wavefunctions for 3 He and 4 He take 
on maximal values if the momenta of the subclusters are zero in their respective rest 
systems. These are for 3 He the momenta of p and d and for He the momenta of 
the two deuterons. This means that for the moving nuclei the subcluster momenta 
should be equal. Therefore to form the a-particle with highest probability in the 
picture of Fig. 1 one has to assume that the two deuterons, d' and d, have equal 
momenta. Likewise for 3 He one has to assume that the proton and deuteron, p and 
d, have equal momenta. This turns out to be kinematically inconsistent. Therefore 
we make a choice and assume that only the two deuterons forming the a particle 
have equal momenta. We justify this choice by the larger binding energy of the a 
particle. 




d 



Fig. 1. Diagram of the reaction mechanism 

It is easy to see that our basic assumption 

k~ d = k d > 



(2-1) 
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fixes the kinematics uniquely. It follows by simple kinematical arguments that 

_, i _, 9 -> -> 

ulab _ ± y cm U lah — h lab (0.0\ 

— 2 P 5 d ~ d V z z ) 

Here the superscripts lab and cm denote the laboratory and 5-body cm systems, 
respectively. Further the total momentum of the picked up proton and the incoming 
deuteron in the lab system is 

K = l% m + pf (2-3) 

Also we get the momentum of the picked up proton in the 3-body center of mass 
system (3CM) as 



kf M = -k™ - -k l d ab (2-4) 



and the 3CM energy as 



■:3CM _ 1 ?cm 3 riab 

Hp 

E 3CM = -!-m CM ) 2 (2-5) 



4m 

We show in Fig. 2 the relevant kinematics for the cm and the 3CM systems. From 
the relation 

kfM = ZpCM _ ^lab (2 . 6) 

it follows under our condition, that the angles shown in Fig. 2 are related as 

q3CM _ q3CM _ q3CM £2.7) 

(note that p = 0p CM = Q c p m ). The dependence of E 3C m on ^ m is illustrated in 
Fig. 3 for 3 deuteron energies. The scattering angle 3CM is shown against 6p m in 
Fig. 4 again for the same 3 deuteron energies. 

c m 





Fig. 2. Scattering angles for the 5-body (cm) and 3-body (3CM) center of mass systems. 

Our claim is now that 0(E d , Q c v m ) w O pd (E 3CM ,9 3CM ) where O pd are the elastic 

pd deuteron polarization observables and O the ones for the reaction 3 He(d, p) 4 He. 

Before calculating these 3N observables we introduce one more approximation. 
Looking at Fig.3 we see that E^cm varies with Op™ and consequently for each 9 p m one 
would have to solve the 3N Faddeev equation. We avoided that for that qualitative 
investigation and have chosen available Faddeev results at three energies which lie 
in the three energy bands for < Q c p m < 40°. They are E 3CM = 66.7, 100, 133 MeV 
corresponding to E d = 140, 200, 270 MeV, respectively. 
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Fig. 3. Effective E^cm energies as a function 
of the proton scattering angle 9p m . The 
solid, dashed and short dashed lines are for 
E d =140, 200 and 270 MeV, respectively. 
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Fig. 4. Effective scattering angle Q 3CM as a 



function of the proton scattering angle 6', 
for the deuteron energies as in Fig. 3 
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§3. Results 



As NN potential we used AV18 in the Faddeev calculations. The operator U for 
elastic scattering has the form (see, for instance, 9 )) 

U = PG 1 +PT (3-1) 

where Go, P and T are the free 3N propagator, permutation operators and a partial 
3N break-up operator, which is determined by a Faddeev equation. The first term, 
the famous nucleon exchange term, is essentially related to the PWIA mentioned in 
introduction. In order to see the importance of solving the Faddeev equation cor- 
rectly and not just replacing U by PGq 1 we compare the corresponding predictions 
for A yy , A xx and A xz in Figs. 5-7. We see large differences especially above about 
15 degrees. Trivially A y is identically zero using only the real term PG^ 1 . 




5 10 15 20 25 30 35 40 



0; m [deg] 

Fig. 5. Tensor analyzing power A yy in elastic 
pd scattering at _Ecm = 133 MeV. The solid 
(dashed) line is calculated from U (PGq 1 ). 
The data point for 3 He(d,p) 4 He reaction is 
from ^ . 
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Fig. 6. The same as Fig. 5 for A xx . 



The predictions of the full Faddeev solution are shown in Figs. 8-11 at £ , 3cm=66.7, 
100 and 133 MeV, respectively. This should be compared to recent data 10 ) . We see 
a behavior qualitatively similar to those data, especially for A y . For the A y data 
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Fig. 7. The same as Fig. 5 for A xz . There is no data point. 
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Fig. 8. The deuteron vector analyzing power 
A y for the pd elastic scattering for 
i?3CM=66.7(solid), 100(long dashed) and 
133(short dashed) MeV as a function of 
9p m , corresponding to the 3 He(d,p) 4 He re- 
action for E d =U0, 200, 270 MeV, respec- 
tively. 
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Fig. 9. The same as in Fig. 8 for A yy . The 
data point for 3 .ffe(d,p) 4 He reaction is from 
i) 
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Fig. 10. The same as in Fig. 8 for A xx . The 
data point for 3 _He(d,p) 4 He reaction is from 



? 0.6 - 




10 15 20 25 30 35 40 
e; m [deg] 

Fig. 11. The same as in Fig. 8 for A xz . 



the minima shift to smaller O^™ 1 value with increasing energy like in Fig. 8. Also 
for Ayy the qualitative behavior is similar in our model and the data, especially at 
the highest energy. For A xx the shapes are again very similar. In Fig. 9 and 10 we 
include one data point from 1 ). This shows that our absolute values are too high. 
For A xz shown in Fig. 11 there are not yet data. 
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§4. Summary and Outlook 

We assumed that the reaction 3 i?e(<i,p) 4 He at forward angles is mainly driven 
by elastic pd scattering. In this model the deuteron picks up a proton from 3 He, 
scatters elastically and combines then again with the spectator nucleons to an a par- 
ticle. Our main assumption is that the momentum of the scattered deuteron equals 
the spectator momentum of the deuteron in 3 He. This leads to a high probability to 
form the final a-particle. The resulting spin-observables are in astonishingly good 
qualitative agreement with the data. Important thereby is, that the elastic pd am- 
plitude is a full solution of the 3N Faddeev equation and not only a simple PWIA 
expression. This model should be generalized by the mechanism that also a neutron 
from 3 He can be picked up. In this case one has to use the nd break-up amplitude. 
Since the polarization of 3 He is carried by more than 90 % by the neutron this second 
mechanism is of course mandatory for a description of C XjX and C y , y . The proton 
pick-up alone is too poor for those spin correlation observables. Also we neglected 
the momentum distributions of the proton in 3 He and of the deuteron in the a par- 
ticle. As an additional improvement the spin of the deuteron should be properly 
rotated for the deuteron polarization observables. 

Based on the promising qualitative results achieved it appears worthwhile to 
improve and enrich the model along the lines mentioned. 
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